Astrophysical constraints on unparticle-inspired models of gravity 
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We use stellar dynamics arguments to constrain the relevant parameters of ungravity inspired 
models. We show that resulting bounds do constrain the parameters of the theory of unparticles, 
as far as its energy scale satisfies the condition Au > 1 TeV and du is close to unity. 
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I. INTRODUCTION 

It has been remarked that the Standard Model (SM) 
is hkely to be incomplete due to the apparent lack of 
scale invariant objects, unparticles [1], besides its well- 
known shortcomings. Implementing scale invariance re- 
quires considering an additional set of fields with a non- 
trivial IR fixed point, the Banks- Zacks (BZ) fields. The 
interaction between SM and BZ fields occurs through the 
exchange of particles with a large mass scale, M*, written 
as 



^BZ = j^OsmObz 



(1) 



where Osm is an operator with mass dimension dsM 
built out of SM fields and Obz is an operator with mass 
dimension cIbz built out of BZ fields. 

At an energy scale Ajj the BZ operators match onto 
unparticles operators (Ou) and Eq. (1) matches onto 
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where du is the scaling dimension of Ou, which can be 
fractional, and Cu is a coefficient function. 

Considering tensor-type unparticle interactions with 
the stress-energy tensor of SM states leads to a modifica- 
tion to the Newtonian potential ^(r), usually referred to 
as ungravity — a gravitational potential with a power- 
law addition [2], 
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where Rq is the characteristic length scale of ungravity, 
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and Au > 1 TeV is the energy scale of the unparticle 
interaction (the lower bound reflects the lack of detection 
of these interactions within the available energy range), 
Mpi is the Planck mass and a is a constant dependent on 
the type of propagator (unity in the case of a graviton). 

The Newtonian potential is recovered for du — ^, 
Rg = (if du > 1) or Rq oo (if du < 1), so that 
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where Rq is the distance where the gravitational potential 
matches the Newtonian one, V{Ro) — $(i?o)- Unfortu- 
nately, the value of Rq is unknown; this may be circum- 
vented by considering only values of du near unity [2] , so 
that Eq. (3) is approximately given by 
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Notice that corrections of this type also arise in the con- 
text of a gravity model with vector- induced spontaneous 
Lorentz symmetry breaking [3] . 



II. POLYTROPIC STELLAR MODEL 

In what follows we examine the bounds on parameters 
Rg and du in Eq. (6) arising from astrophysical consid- 
erations about stellar equilibrium. In order to do so, we 
shall extend considerably the range of ungravity correc- 
tions. Before discussing these bounds in detail we point 
out that astrophysical and cosmological constraints on 
unparticles have been discussed in Refs. [4-9], and the 
ones arising from nucleosynthesis have been studied in 
Ref. [10]. We also mention that the technique to be 
employed has been developed to constrain Yukawa type 
corrections to the Newtonian potential [11] as well as 
to examine alternative gravity models with nonminimal 
coupling between curvature and matter [12]. 

The simplest model available for stellar structure in- 
volves the polytropic gas model, which assumes the state 
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equation P = _fsrp("+^)/", where P is the pressure, p is the 
density, n is the so-called polytropic index and K is the 
polytropic constant. The above equation of state allows 
one to write the relevant thermodynamical quantities as 

p = Per(0 , T = TM) , P = Pcr+i(0 . (7) 

where pc, Tc and Pc correspond to the values of density, 
temperature and pressure at the core of the star, respec- 
tively. The dimensionless function 9{£,) depends on the 
dimensionless variable ^, related to the radial coordinate 
through r = (3^, where 
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Using Eqs. (7), the hydrostatic equilibrium condition 

-^f^fU-G^ (9) 
ar \ p ar J ar 



may be rewritten as 
1 d 



(10) 



the Lane-Emden equation [13]. This differential equation 
is subjected to the initial conditions 9{0) = 1 and 9'{0) = 
0. A solution to the Lane-Emden equation allows for 
the determination of the thermodynamical quantities of 
a star in terms of their values at its center. The profile of 
6{^) depends only on the choice of the polytropic index 
n, not on the size of the star, manifesting the homology 
symmetry of this equation. 



III. MODIFIED LANE-EMDEN EQUATION 

In this section, we develop a method similar to that 
presented in Rcf. [11], in order to extract the relevant 
bounds on du and Rq- We consider the modified poten- 
tial Eq. (6) and assume the validity of the Newtonian 
regime (low density and small velocities) to obtain the 
modified hydrostatic equilibrium equation: 
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After some algebraic manipulation, this can be cast as 
1 d (r'^dP(r) 
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Performing the substitutions r = and p = PcO^, we 
obtain the perturbed Lane-Emden equation 
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where = Rg/I^ has been defined, for convenience. 
Using relation dM{r)/dr = 4TTp{r)r'^, together with def- 
initions Eqs. (7) and (8), we obtain 
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(14) 



which can be used to simplify the second term on the 
r.h.s. of Eq. (13), which now reads 
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It is interesting to point out that this modified Lane- 
Emden equation, unlike Eq. (10), has no homology sym- 
metry, due to the presence of in Eq. (15) — and hence 
the stability of the star will depend on its radius. The 
unperturbed central temperature Tc of a star is obtained 
from the solution 9o{0 of Eq. (10) [13], 
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where signals the surface of the star, defined as 
^o(Cio) = 0- Considering n = 3, which describes fairly 
well the overall features of the Sun, one finds ~ 6.90 
[13]. In the presence of the ungravity perturbation into 
the gravitational potential Eq. (6), the central temper- 
ature Tc will be shifted from Tco, the value obtained by 
using the solution 9o to the unperturbed Lane-Emden 
equation, Eq. (10), yielding the ratio 



(17) 



We now seek a numerical solution of Eq. (15) that allows 

us to estimate the ratio Eq. (17) to obtain a contour plot 
of the relative shift T^ — 1 for different values of djj and 
Rg (for n — 3). 

We consider two ranges of values for du and Rg- du ^ 
1 for Rg < Rs and du < 1 for Rg > Rs, Rs 10^ m 
being the Sun's radius. For du si 1, we also assume 
that the length scale Rg is larger than the Schwarzschild 
radius of the Sun, Rm = 2GMq/c^ k. 1.5 km, so that 
no relativistic corrections of the form Rm / r have to be 
considered. In what concerns the modified Lane-Emden 
equation, Eq. (15), the following ranges are considered: 
1.0 < du < 1.06 for < < 1 and 0.94 < du < I for 
10 < < 10^. 



IV. RESULTS 
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FIG. 1: Contour plot of T,. — 1 in function of Ra and du- 




FIG. 3: Lower bound of log(M,/Mpi) for a = (black curve), 
a = 2/3 (dark grey curve), a = 1 (light grey curve), a = 1.9 
(black, lower curve), and Au = 1 TeV (solid curve), Au = 
10^ TeV (dashed curve). 
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FIG. 2: Contour plot of — 1 in function of log Ra and du- 

Numerical solutions of Eq. (15) allow for obtaining 
contour plots for T,. — 1 as a function of du and ^g- The 
results are depicted in Figs. 1 and 2 for jT^ — 1| < 0.06, 
the uncertainty in the Sun's central temperature [13]. 

Designating the line in Fig. 1 that indicates a 6% 
change of the Sun's central temperature as i?_(d[/), one 
sees that Rcidu) > R-{du) for du < 1. Thus, Eq. (4) 
leads to 



^ > [nAuR-idu)]'-'''' f{du,a) 
Mpi 



where 
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is defined, for convenience. One might plot the lower 
bound obtained above as a function of du, fixing the 
model parameters a and Au- This is depicted in Fig. 3, 
for a = 0, 2/3, 1, 1.9 and suitable values for Au; all 
lines converge to the trivial point du = I, > 0. 



Similarly, one obtains from Fig. 2 the upper bound 

Rcidu) < R+{du), where the latter denotes the line cor- 
responding to the 6% change in the Sun's central tem- 
perature. Resorting again to Eq. (4), this again yields 

^* >[nAuR+{du)]'-''" f{du,a) . (20) 



log(M./Mp,) 



Mpi 



1.01 1.02 1.03 1.04 1.05 1.06 



du 




FIG. 4: Lower bound of log(M*/Mpi) for a = (black curve), 
a — 2/3 (dark grey curve), a = 1 (light grey curve), a = 1.9 
(black, lower curve), and Au = 1 TeV (solid curve), Au = 
10^ TeV (dashed curve). 

The obtained lower bound is depicted in Fig. 4 (as 
before, the lines converge to the point du = I, M* > 0). 



A. Discussion 

As stated in Ref. [2], corrections to the Newtonian 

potential $Ar(r) with du < ^ might appear unfeasible, 
since these will overcome the 1/r dependence of $Ar(r) 
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for r > Rg and lead to long-range deviations. This might 
be alleviated by letting Rq be so large that this crossover 
occurs well beyond the relevant astrophysical range, and 
one can no longer assume a static, spherically symmetric 
Ansatz for the metric. 

Alternatively, one may consider values so close to unity, 
du ^ ^1 that the perturbation to the Newtonian po- 
tential, Eq. (6) may bo expanded as {r/RcY~'^'^" « 
1 + 2{du — 1) \og{r / Rg) , and the logarithmic dependence 
is attenuated by the small djj — 1 term: for instance, 
for r = 100 AU, the typical dimension of the Solar 
System and Rq ~ 10 AU, the maximum value consid- 
ered here, this yields {r/Rcf''^'^" k. I + 4{du - 1); 
assuming the same value for Rg and instead, if r ~ 
100 kpc ~ 10^° AU the size of a galaxy, one still gets 

(r/i?G)'-"^^ «l + 20(dc/-l). 

With these considerations in mind, the method de- 
veloped here shows that one can successfully constrain 
the range of Mu for rfj/ < 1: in particular, assuming 
Ajj > 1 TeV one achieves lower bounds ranging from 
> (10~^ — 10^)Mp( (even lower bounds can be ob- 
tained for values of du closer to unity). 

For the case djj ^ 1, one obtains a lower bound ex- 
hibiting a peak around djj ~ 1.01, with typical values 
Mh, > (10^2 _ lO"i)Mpi. Rcf. [2] presents lower boimds 
for as a function of Au, for d;/ = 2, 3, 4 — values 
which are beyond the reach of this study. In a subsequent 
study, the cases du = 1.1, 1.5, 2 were considered, with 
the first case closer to the range considered here [9] . 

By solving Eq. (15) for rfc/ = 1.1 and finding the value 
of Rg that yields T^ — 1 = 0.06, one obtains a lower bound 
of about > (lO-'^ - 10-2)Mp( w (lO^^ - lO^") TeV, 
depending on Au and M* (this may be checked by ex- 
trapolating the plot of Fig. 2). This limit is much 
greater than the one foimd in Rcf. [9], where a result 
Af, > (10^ - 10^) TeV is reported (for Au = 10** TeV). 
This indicates that the developed method hints at a much 
more stringent bound for M*, for d[/ > 1. 



static equilibrium conditions within a star. This leads to 
a perturbed Lanc-Emden problem that is then examined 
for the polytropic index n = 3. Prom the resulting change 
in the star's central temperature, we obtain constraints 
on the ungravity parameters Rg and du- Given that 
the overall properties of the Sun are well described by 
the n = 3 case, we allow for the ungravity correction to 
affect this up to the upper bound on the Sun's central 
temperature, ATc/Tc » 0.06. 

We find that, for du > I and Au > 1 TeV, lower 
bounds on M, are in the range (10^^ — 10^^)Afp;. For 
du ^ ^ and Au > 1 TeV, must lie in the range 
above (10~^ — 10^)Mpi. Of course, our bounds are com- 
plementary to the ones obtained from torsion balance 
experiments which test a much smaller range of Rq [14] , 
actually about 80 /xm. 

The reported results for > 1 are either more strin- 
gent [4, 5, 7-9] or similar [6] to those previously available. 
The lower bound derived for du ^ 1 is more relevant, 
since it has been not obtained so far. In our opinion, 
this arises from misconception that a negative exponent 
in Eq. (6) is disallowed by long-range experiments [14]: 
while this is true for large values otl — du, a range closer 
to unity, dt/ < 1 yields an approximately logarithmic 
correction, with large deviations from the Newtonian po- 
tential suppressed by the smallness of the c?c/ — 1 term. 
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